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Abstract 



We present a mean field treatment of a strongly correlated model of elec- 
trons in a three-dimensional anisotropic system. The mass of the bare elec- 
trons is larger in one spatial direction (the c-axis direction) , than in the other 
two (the a6-planes). We use a slave fermion decomposition of the electronic 



Q ■ degrees of freedom and show that there is a transition from a deconfined to a 

^ ■ confined phase in which there is no coherent band formation along the c-axis. 



One of the most controversial, and hard to understand, problems related with high-T^ 
cuprates is the anomalous charge transport observed experiment allytl. The charge dynamics 
^ I reflects the anisotropy in the crystal structure of these compounds, which consists of weakly 

O ' coupled planes. In the usual notation, we will refer to "c-axis" and "a6-planes" as the 

^ ■ directions transverse and parallel to the planes respectively. The in-plane conductivity aab 

■^ ! shows a behavior characteristic of the metallic state. On the other hand, close to the insu- 

J3 \ lating state, in the so called underdoped regime, the c-axis conductivity a^ is "incoherent" : 

the values of ac are below the minimum metallic conductivityH, the temperature depen- 
n^ • dence is anomalous, and the frequency dependence does not show signatures of Drude-like 

^ ■ behaviorH'Q. 

Q , Band structure calculations indicate an anisotropy which, within the framework of Boltz- 

man transport, imply metallic behavior with an anisotropy cXc/cTab well above the experimen- 
tal observation. Perturbative treatments within the Fermi-liquid theory indicate that the 
anisotropy is not renormalized by interactions^. Perhaps the main objection to the "conven- 
tional" theories of c-axis transportQ'B is the observed value of the anisotropy of the conduc- 
tivity. In the superconducting phase coherence is reestablished in all directions^. This lead 
Anderson and others to attribute the anomalies in transport in the normal state to the effect 
of strong electronic correlations, and to conclude that in order to describe the incoherent 
c-axis conductivity the Fermi-liquid picture should be abandoned. The starting point used 
as a paradigm is the one-dimensional correlated problem where it is rigorously known that 
the Fermi-liquid picture fails. Considerable work has been done in weakly coupled chains 
that suggest that a state can be formed in which the coherence is confined to the motion 
along the chains, the motion transverse to the chains being incoherentH. 



A complete theory for the charge dynamics in anisotropic strongly correlated systems is 
not yet available. Due to the complexity of the problem, much work remains to be done 
in order to develop a fully consistent and controlled calculation scheme that could account 
for the phenomenology indicated by the experiments. In the mean time, the analisys of 
simple models is useful as a starting point towards the final answer. Here we present a 
mean field treatment of a system of coupled planes that includes the strong anisotropy, and 
incorporates the strong correlations responsible for the non-Fermi liquid behavior. We show 
that, within that mean field, a transition from a deconfined to a confined phase takes place. 
The parameter signaling the transition is the gain in kinetic energy due to band formation 
in the c-axis direction. 

We consider the Hubbard model in the limit of infinite on-site repulsion, described by 
the following Hamiltonian 

H =Y1 kj 11(1 - ni,_^)4^Cj- ^(1 - n^- _^), (1) 

where (i, j) refers to near neigbors on a cubic lattice where the anisotropy is incorporated 
in the values of hopping the matrix elements: tj j = t\\ for in-plane hoppings and tjj- = t^ 
for the motion along the c-axis. The fermion operators c\„ create an electron at site i only 
if the site is empty. 

A well known mean field description of Hamiltonian |l| is the slave-bosonllylld approach 
in which each local configuration has associated with it a fermionic or bosonic degree of 
freedom, such that q „. = alg-Cj, where a\^^ creates a fermion with spin a at the z-th site 
representing a singly occupied configuration, and Cj destroys a boson representing the empty 
state at the same site. A standard mean field calculation decouples fermions and bosons 
and relaxes the exact constraint of one "particle" (fermion + boson) per site. The resulting 
problem is that of non-interacting bosons and fermions self-consistently coupled. As a result 
the ideal bosons condensate in a A; = state, the overall effect being a renormalization of 
the masses of the fermions. It is important to note that, even for an anisotropic system, 
the fc = bosonic ground state wave function does not know about the anisotropy, and the 
mass renormalization is the same in all spatial directions. Consequently, such an approach 
preserves the anisotropy and the Fermi liquid character of the ground state. At least formally, 
one can conceive corrections to this state that improve the treatment of the constraint to 
avoid multi-occupancy of the particles at the same site. There are, however, other alternative 
treatments that-still within the mean field level-take into account the hard core constraint 
for the bosons exactly. In the present work we present a mean field along this line. In 
what follows we show that at an alternative description in terms of slave-fermions for the 
infinite-f/ case breaks the Fermi liquid description and produces a confined coherent satellHl 
in the ah planes. 

We introduce a description in which the original projected fermions are represented by 
three fermions: 



= C,; 



-'^,c^ — ^i.a 



r{l - Ui^-^) = tti^^fl^fi^l (2) 



The above representation respects the anticommutation relation between the projected 
operators Ci^a and c|^ provided one stays within the physical Hilbert space. A related fermion 
linearization was presented in Ref.li^. 



The product fi^^fi,i is a spin flip operator corresponding to a pseudo-spin degree of 
freedom not related to a. When this fictitious spin is J, in site i this means that the site is 
occupied, and the site is empty if the spin is |: there are as many /j^'s as there are electrons 
and as many /^'s as there are holes. The / fermions therefore satisfy 

(4tAt) + (fliAi) = 1, E(«.W) + (fliki) = 1, (3) 

and, in turn 

Y,{4,a(^i,a) = I - S, (4) 

(T 

with S representing the fractional deviation in occupation number with respect to the half- 
filling case of one electron per site. 

At the mean field level the ground state wave function consists of a direct product of 
three Fermi seas, one per each of the fermion degrees of freedom. The total energy in this 
approximation is given by 

Eo = -T.kjA3xlp (5) 

(id) 

with 

Aij = Y,{4,a(^j,^), (6) 

cr 

X., = (fy,,) = (fiki)^ (7) 

where the last equality holds because we are dealing with a bipartite lattice with particle-hole 
symmetry. The three species of fermions are free with their hopping amplitudes renormalized 
by the factors Aij, and Xi,j- These factors are responsible for renormalizing the anisotropy 
and can be better visualized in the mean field Hamiltonian 

HmF = — 2_/ ^«J Z2 [Xi,jO'i,CTO'j,(7 + A.ijXi,jfi,CTfj,(7 + C, (8) 

with C a constant. 

Note that, for small deviations from half filling, the f\^ {fi,i) fermions are moving close 
to the bottom (top) of their band, whereas the a fermions are close to the center of the 
band. This makes their respective Fermi surfaces different. 

Our mean field can be understood in two steps. First the a fermions are decoupled 
from the / fermions. At that level, the a fermions are free, but the /| fermions and the fi 
fermions are strongly correlated. The dynamics of the system of / fermions at this level is 
identical to that of an xy-model, and can be mapped onto a hard-core boson problem. In 
the second step the / fermions of different spin are decoupled and treated as free fermions 
(with a self consistent constraint on the dynamics)lij. Note that, at the level of step one, 
the above mentioned system of hard-core bosons will in principle have an anisotropy in the 



expectation values of the kinetic energy terms that will depend on direction. This is due to 
the quantum fluctuations introduced by the hard-core constraint. 

At half fllling (6 = 0), the kinetic energy of the / fermions is zero, the renormalization 
factor Xij = 0, giving the localized limit of the a fermions which we identify as the Mott 
insulating state. On the other hand, far from half fllling, for 5 ~ 1, the density of a -fermions 
is so low that they should behave as non-interacting, but our mean fleld fails to recover this 
limit. Decoupling the /-fermions from the a-fermions is not a good approximation in the 
limit of high doping S because the probability of flnding an /^-fermion at a site occupied 
by an a fermion is very low [~ (1 — S^], while the exact dynamics requires this probability 
to be one. Therefore our results will be valid close to half fllling, or 5 ~ 0. 

Due to translational invariance, the ground state energy will be a function of the four 
quantities A\\, A±, x\\ and x±'- 

Eo = -4t||A||xJ - 2t^A^xl- (9) 

The one particle energies of the / and a fermions are respectively 

eI = t||A||X||£fe|| -t±Aj_x±'2coskz, (10) 

^k = ^||x5^fc||-t^xl2cosA;„ (11) 

with 

Ekn = —^{cos kk + cos ky). (12) 

Effective chemical potentials A and fi have to be determined for each of the two types of 
fermions through the equations 

^E/(^i^-^) = '^' ^E/(^k-/^) = ^^- (13) 

We approximate the reduced density of states corresponding to the motion within the 
plane by a constant: J2ku ^{^ ~ ^fe,,) = ©(4 — |£|)/4, and flnd that the mean fleld equations 
can be written in terms of the parameters a and (3 deflned as 

After straightforward integrations, and using the fact that close to half fllling the Fermi 
surface of the a fermions is open, the mean fleld equations are 

1.. ,2^ ff^V . P 
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-4,| = 2(1-*V i . ^x = ^. (15) 



2kQ- — \kQ + -sm.2kQ\ - —sinko) , (16) 




X± = ^ I (l - ^J 2/cosinfco + f (^0 + ^ sin2A;o) 1 , (17) 

with A = A/(t||yl||X||) determined from the equation 

S = — {{X + 4:)h + 2asmko}, (18) 

and ko = cos"^[— (A + 4)/2a] for |(A + 4)/2a| < 1 and vr otherwise. Note that a plays the role 
of an effective anisotropy of the /-fermions. For a given 6, if we fix a, the renormalization 
factors X ci^nd A are determined by the Equations (|T5|) through (|1^). This means that a 
plays the role of a variational parameter with respect to which we have to minimize the 
energy Eq. As an example, in Figure ^ we show some curves of Eq vs. a for different values 
of doping using as a parameter the bare anisotropy t±/t\\. 

The curves indicate that for fixed t±/t\\ there is a discontinuous jump in the position of 
the minimum of Eq as 6 is varied. The curve shown in Figure |^ for 6 = 0.002 corresponds 
to the confined phase for which a = 0, and the renormalization factor x± = 0. The curve 
for 6 = 0.0018 has its minumum at finite a and hence corresponds to a three dimensional 
metal with a renormalized anisotropy. 

A phase diagram that result from our calculation is shown in Figure ^ 

A very important point is to establish that the particle motion does not correspond to a 
Fermi liquid. We show this by computing the form of the occupation number of the original 
fermions in the confined phase within our mean-field squeme: 

n.,. ^ (etc.) = ^ + ^Ee''^'^'-^^^(4c,,.) 

with n the particle density. The term {cj^Cj^^) is evaluated using the representation of 
Equation |^. In mean field the result is a convolution of the occupation numbers of the three 
fermions (/'s and a). Using the constraints of Equations ^ and ^ one obtains 

^^ pq 

The occupation numbers above correspond to three Fermi surfaces. For small 5 the 
Fermi surfaces corresponding to the / fermions are are two circles centered respectively at 
k = (0, 0) and k = (vr, vr). On the other hand, the Fermi surface of the a fermions are close 
to a diamond. The result of the convolution above (See Figure 4) is that Uk^a does not have 
a discontinuity, implying an non-Fermi liquid state. 

A few points related to the calculation deserve a comment: 

i) Due to the approximation made in the density of states we cannot recover the isotropic 
case. The approximation used is aimed at describing anisotropic systems. 

a) In our calculation the confined regime is identified by the vanishing of the expectation 
value of the interplane hopping indicating that there is not band formation along this direc- 
tion. We interpret this result as an indication of incoherence, even though one expects some 
interplane-coupling to remain in the exact incoherent regime. The picture is analog to the 



slave boson description of the Mott insulator. There, the insulating state is characterized 
by a vanishing of the inter-site hopping, while we know that in the exact ground state this 
magnitude is small but finite. 

In summary, we have presented a mean field calculation and derived a phase-diagram of 
a strongly interacting anisotropic system. We have shown that, as the anisotropy increases, 
for small deviations from half filling a transition takes place from a deconfined phase to a 
confined phase in which the motion in the c-axis direction is completely incoherent while 
the motion in the afe-direction corresponds to a coherent, non-Fermi liquid state. 
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FIG. 1. Schematic rendition of the slave fermion decomposition. When an a fermion hops 

from site i + 1 to site i, there is a spin-flip of / fermions represented by the dashed arrows. The 

upper (lower) part of the figure represents the configuration before (after) the hopping process. 
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FIG. 2. Ground state energy vs. the variational parameter afor an anisotropy of tx/tii = 0.3 
and the indicated values of doping. 
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FIG. 3. Phase diagram valid in the low doping regime indicating the boundary between a 
confined phase and a three-dimensional anisotropic phase. 
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FIG. 4. Fermion occupation number vs. wave vector in the direction marked in the inset. The 
dashed hne shows the bare, non-interacting, occupation number. Note that there is no discontinuity 
in nfc 0-) as expected in a non~Fermi-hquid state. The inset also shows the Fermi surfaces-in the 
first quadrant only-of the / fermions (smaU circle shown in short dashed and arc in continuous 
line close to the M point) and of the a fermions (dashed line). 
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